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Abstract 

In this paper, we are interested in the numerical solutions of stochastic functional 
differential equations (SFDEs) with jumps. Under the global Lipschitz condition, we 
show that the pth moment convergence of the Euler-Maruyama (EM) numerical solu- 
tions to SFDEs with jumps has order 1/p for any p > 2. This is significantly different 
from the case of SFDEs without jumps where the order is 1/2 for any p > 2. It is there- 
fore best to use the mean-square convergence for SFDEs with jumps. Consequently, 
under the local Lipschitz condition, we reveal that the order of the mean-square con- 
vergence is close to 1/2, provided that the local Lipschitz constants, valid on balls of 
radius j, do not grow faster than logj. 

Keywords: Euler-Maruyama; Local Lipschitz condition; Stochastic functional differ- 
ential equation; Rate of convergence; Jump processes. 
Mathematics Subject Classification (2000) 65C30, 65L20, 60H10. 



1 Introduction 

Recently, the theory of functional differential equations (FDEs) has received a great deal 
of attention. Hale and Lune ^ have studied deterministic functional differential equations 
(DFDEs) and their stability. For stochastic functional differential equations (SFDEs), we 
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here highlight the great contribution of Kolmanovskii and Nosov [7j and Mao [S]. Kol- 
manovskii and Nosov [7| not only established the theory of existence and uniqueness of 
SFDEs but also investigated the stability and asymptotic stability of the equations, while 
Mao [8] studied the exponential stability of the equations. 

On the other hand, stochastic differential equations (SDEs) with jumps have been widely 
used in many branches of science and industry, in particular, in economics, finance and 
engineering (see, for example, Gukhal [3], R. Cont Sobczyk [TT] and references therein). 
Since most SDEs with jumps cannot be solved explicitly, numerical methods have become 
essential. Under the local Lipschitz condition, Higham and Kloeden |3] showed the strong 
convergence and nonlinear stability for the EM numerical solutions to SDEs with jumps, 
while, in [5], Higham and Kloeden further revealed the strong convergence rate for the 
backward Euler (BE) on SDEs with jumps, provided that the drift coefficient obeys one-side 
Lipschitz condition and polynomial growth condition. 

Returning to the SFDEs, we recalled Mao [9] developed a numerical scheme for them. 
Under the local Lipschitz condition, Mao [9] showed the strong convergence of the EM numer- 
ical solutions to SFDEs, but revealed the rate of the convergence under the global Lipschitz 
condition. But there is so far no work on numerical methods for SFDEs with jumps. 

Motivated by the papers mentioned above, we are here interested in the numerical so- 
lutions to SFDEs with jumps. In comparison with the results obtained by Mao [9], our 
significant contributions are: 

• Under the global Lipschitz condition, we show that the pth moment convergence of 
the EM numerical solutions to SFDEs with jumps has order 1/p for any p > 2. This 
is significantly different from the case of SFDEs without jumps where the order is 1/2 
for any p > 2. In practice, it is therefore best to use the mean-square convergence for 
SFDEs with jumps. 

• Under the local Lipschitz condition, Mao [9] showed the strong convergence without 
rate of the EM numerical solutions to SFDEs without jumps. However, we shall reveal 
that the order of the mean-square convergence is closed to 1/2, provided that the local 
Lipschitz constants, valid on balls of radius j, do not grow faster than logj. More 
precisely, the order of the mean-square convergence is 1/(2 + e), provided that the 
local Lipschitz constants do not grow faster than {logjY^^^~^^\ 

• Some new techniques are developed to cope with the difficulty due to the jumps. 

This paper is organized as follows: Section 2 gives some preliminary results, in particular, 
the EM numerical solutions to SFDEs with jumps are set up. In section 3, we discuss the 
pth moment convergence of the EM numerical solutions to SFDEs with jumps under the 
global Lipschitz condition. The rate of the mean-square convergence of the EM numerical 
solutions to SFDEs with jumps under the local Lipschitz condition is provided in Section 
4. Finally, in order to make the paper self-contained, an existence-and-uniqueness result of 
solutions to SFDEs with jumps is provided in the Appendix. 
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2 Preliminaries 



Throughout this paper, we let {Q, J-", {J-'t}t>o, P} be a complete probability space with 
a filtration {J-'t}t>o satisfying the usual conditions (i.e., it is continuous on the right and 
jFo-contains all P-zero sets). Let | ■ | denote the Euclidean norm and the matrix trace 
norm. Let r > and D := D{[—t, 0]; i?") denote the family of all right-continuous functions 
with left-hand limits ip from [— t, 0] to i?", and D := £)([— r, 0]; i?") denote the family of 
all left-continuous functions with right-hand limits from [— r, 0] to i?", we will always 
use llv^ll = sup_^<5i<Q \(p{9)\ to denote the norm in D and D potentially involved when no 
confusion possibly arises. D'^^^{[—t,0];R'^) denotes the family of all almost surely bounded, 
jFg-measurable, D([— r, 0]; -R"')-valued random variables. For all t > 0, Xt = {x{t + 6) : 
—T < 6* < 0} is regarded as a D([— r, 0]; i?"')-valued stochastic process. Let x(t~) denotes 
lims|tx(s) on t G [— t, T] and Xt- = {x(t + 9)~ : — r < 9 < 0}, it is easy to see x{t~) is a 
D{[—T, 0]; -R"')-valued stochastic process. 

In this paper, we consider the following SFDE with jumps 

dx{t) = f{xt-)dt + g{xt-)dB{t) + h{xt-)dN{t), 0<t<T, (2.1) 

with the initial data xq = ^ E D^^([-r, 0]; P"). Here, f,h : I)([-r, 0]; i?") ^ R"", g : 
D[[—r,0]]R"') — > i?"^™, B{t) is an m-dimensional Brownian motion and N{t) is a scalar 
Poisson process with intensity A. We further assume that B{t) and N{t) are independent. It 
should be pointed out that the solution of Eq. (12. ip is in D{[—t, 0]; i?") 

For our purposes, we need the following assumptions which can also guarantee the exis- 
tence and uniqueness of solution to (12. ip (see Appendix). 

(HI) (Global Lipschitz condition) There exists a left-continuous nondecreasing function fi : 
[-r, 0] ^ R+ such that, for all e D{[-t, 0]; P"), 

\f{^) - fm' V \g{^) - gm' V iM^) - < J' - mi'dm- (2.2) 

Remark 2.1. For simplicity, we write L = /i(0) — ix{—t), which is referred to as the global 
Lipschitz constant. Note from (12. 2p that, for all 1^,1^ E D{[—t, 0]; R^), 

|/(^) - fm' V \g{^) - gin' V \h{^) - <Ly- • (2.3) 

This further implies the linear growth condition; that is, for ip G /)([— r, 0]; -R"), 

|/(y.)pV|^7(y.)pV|M^)P<K(l + ||v.f), (2.4) 

where K = 2{Ly |/(0)p V |^(0)|2 V \h{Q)\'^). 

(H2) (Continuity of initial data) For ^ G D^^^{\-t, 0]; P"), < m < r and j9 > 2, there is a 
constant /3 > such that 

E\ sup |e(.)-e(t)r I (2.5) 

\ -T<3<t<0. / 
\ \t-3\<U I 
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For given T > and r > 0, the time-step size A G (0, 1) is defined by 

N M 

with some integers N > t and M > T. The EM method apphed to ( 12. ip produces approxi- 
mations y{kA) ^ x{kA) by setting y{kA) = ^(fcA), — < k < 0, and 

^((A; + 1)A) = y{kA) + fiyuA) A +g{yki.) ABk + h{ykA) A Nk, (2.6) 

where ABk = B{{k + l)A)-B{kA) is a Brownian increment, ANk = N{{k + l)A)-N{kA) 
is a Poisson increment, and ytA = {ykAid) '■ —t < ^ < 0} is a D{[—t, 0]; i?")-valued random 
variable defined by 

VkAie) = il±il^2/((fc + ^) A) + ^-^y{{k + I + 1) A) (2.7) 

for iA < 9 < (i + 1)A, i = —N, —{N — 1), ■ ■ ■ , —1, where in order for |/_a to be well defined, 
we set yi-iN + 1)A) = ^(-iVA). 

Given the discrete-time approximation {y{kA)}k>Q, we define a continuous-time approx- 
imation y(t) by y(t) = ^(t) for — r < t < 0, while for t G [0,T], 



yit)=m+ ff{ys-)ds+ f g{ys-)dB{s) + [ 
Jo Jo Jo 

where, for fixed 9 G [— r, 0], 

M-l 

yt- = limys, Vt = l/fcA/[fcA,(fc+i)A)(^)- 



\{y,-)dN{s), (2. 



k=0 

It is easy to see y{kA) = y{kA) for k = —N, —N + 1, ■ ■ ■ , M. That is, the discrete-time and 
continuous-time EM numerical solutions coincide at the gridpoints. 

Remark 2.2. It is easy to observe from (12. 7p that 

\\ykA\\= max \y{{k + i)A)\, k = -1,0,1, M - 1, (2.9) 

—N<i<0 

which further yields 

WukAW < \\ykA\\, k = -1,0,1, - ■ ■ ,M - 1, 

by y{kA) = y{kA) and, for any t G [0,T], 

ll^tll = \\y[jr]A\\ < \\y[i]A\\ < sup \y{s)\, (2.10) 

~T<S<t 



where [^] is the integer part of ^. 
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3 Convergence under the global Lipschitz condition 



In this section, we will investigate the rate of the convergence under the global Lipschitz 
condition. Our results reveal a significant difference from these on the SDEs without jumps. 

Lemma 3.1. Under the condition l[2.4\) , for any p > 2 there exists a positive constant 
H{p) := H{p, T,^,K) which may dependent on p,T,^, K such that 

e( sup \x{t)\p)yE( sup \y{t)A < H{p). (3.1) 

V -T<t<T ' \-T<t<T J 



Proof. Since the arguments of the moment bounds for the exact and continuous approx- 
imate solutions to (12.11) are very similar, we here only give an estimate for the continuous 
approximate solution y{t). For every integer R> 1, define the stopping time 

Or = mi{t > : \\yt\\ > R} 

It is easy to see from (12. 8p that, for any t E [0, T], 



E[ sup \yisAeR)-\n <E{ snp \yisAeR)f 

yO<S<t I \0<S<t 



EUr + E{ sup 

0<s<t 



f{y{rAeR)-)dr 



(3.2) 



+ E\ sup 

vO<s<t 



9{y{rAeR)')dB{r) 



+ E ( sup 

,0<s<t 



Noting that E\\y(^if^gj^y\\ < R and (I2.10p . one may have E\\y(^if^g^y\\ < R. By the Holder 
inequality and (12.41) . 



E \ sup 

^0<s<t 



f{y{rAeRr)dr 



<TP-' f E\f{y^rMn)-)\'dr 
Jo 



'pdr 



(rAdnrW^Pdr 



<TP-i / E[K{l + \\y^r 
Jo 

= T^-' [ E[K{l + \\y^, 
Jo 

^0 

This, together with (I2.10p . immediately reveals that 



E ( sup 

yO<S<t 



f{y(rAen)-)dr 



<ciT + ci E[ sup \y{r A9Ry\P ] ds, 

Jo \—T<r<s 



(3.3) 



where Ci = 22~'Tp-^K2. Now, using the Burkholder-Davis-Gundy inequality [8, Theorem 
7.3, p40] and the Holder inequality, we deduce that there exists a positive constant Cp such 
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that 



E [ sup 

.0<s<t 



g{y(rAen)-)dB{r) 



P\ / rt \P/2 

<c„E i I \g{y{rAen)-)\'^dr ) 







< CpT'-^ / E\giy. 
Jo 

In the same way as (13.31) was done, it then follows easily that 



\Pdr. 



E ( sup 

vO<s<t 



9{y{rAen)-)dB{r) 



]<C2T + C2 [ e( sup \y{rAeR)-Ads, 

J Jo \-T<r<s / 



where C2 = K^Cp. Moreover, observing that N{t) = N{t) — At,t > is a martin- 

gale measure, using the Burkholder-Davis-Gundy inequality fUlj Theorem 48, pl93]. Holder 
inequality and (12 ■41) . we obtain for some positive constant Cp, 



E ( sup 

^0<s<t 







<E( sup 

^0<s<t 



h{y{rAeR)-)dN{r) 



Hy{rAen)-)dN{r) + X / h{y (^rAe a)-) dr 
^0 



< 2P 

< 2^ 



E [ sup 

0<s<t 



h{y(^rAeR)-)dN{r) 



+ \^ sup 

0<s<t 



h{y{rAeR)-)dr 



p/2 



where = 2^2 ^K^ 



c,A^/^S^ \h{y(rAeRy))?dr^ + X^T^-' / E\h{y(^r^eR)-)Wdr 
<c^T + c^f E^sup \y{rAeRYAds, 

Jo \-T<r<s / 

. Hence, in 



\0<s<t 

< 4P-^ 



E[ sup \y{sAeR)-\P 

EUr + {ci + C2 + C3)T + {c, + C2 + C3) [ e( sup lyirAOt 

Jo \-T<r<s 



ds 



Note that 



E[ sup \yisAeR)-\P\ <EU\\P + E{ sup IvisAOR) 

^T<S<t J \0<S<t 



Applying the Gronwall inequality and letting i? ^ 00, we then obtain 



Ei sup \yin\P\ <H{p). 

\-T<t<T / 

Since T is any fixed positive number, the required assertion follows. 

In order to obtain our main results, we need to estimate the pth moment of y{s- 



□ 
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Lemma 3.2. Let the conditions and \2. 5]) hold. Then, for p > 2 and s G [0,T], 



+ <7A, -r<e<0, (3.4) 

where 'j is a positive constant which is independent of A . 

Proof Fix s e [0, T] and 9 e [-r, 0]. Let ks G {0, 1, 2, ■ ■ ■ , M - 1}, ke E {-N, -N + 
1, ■ ■ ■ , —1} be the integers for which s G [kgA, {ks + 1)A), 9 G [kg A, {kg + 1)A), respectively. 
For convenience, we write v = s + 9 and k^ = kg + kg. Clearly, < s — kgA < A and 
< ^ - fceA < A, so 

0<v- k^A < 2 A . 
Recalling the definition of ys, s G [0,T], we then yield from (12. 7p that 

ysi9) = yk.Ai9) = y{k.,A) + ^'^^ [y{{K + 1)A) - y{KA)l 

which implies 

E\y{s + 9)- y,(e)r < 2^-^E\y{{k, + 1)A) - y{k,A)\P + T'-^E\y{v) - y{k,A)\P . (3.5) 

For k^ < —1, it thus follows from (12.51) that 

E\y{{k, + 1)A) - y{k,A)\^ < (3 A . (3.6) 

Note that, for some H := H{m,p), 

E\B{t)\P < Ht^, t>0 (3.7) 

and, by the characteristic functions argument, for A G (0, 1), 

E\ANk\''<CA, (3.8) 

where C is a positive constant which is independent of A. For ky > 0, using (12. 6p and noting 
giVkyA) and Bk^, h^yk^^) and N^^ are independent, respectively, we compute 

E\yiikv + l)A)-yik,A)\P 

< 3^-' [E\fiy,^^)\P A^ +E\giy,^^)\PE\ A + E\hiy,^^)\PE\ A N.J^] . 

Taking (12.41) into consideration and applying Lemma [3.11 we then obtain that for A G (0, 1) 

E\y{{k^ + 1)A) - yik,A)f < 3P-^2-2-^K^l + H{p)){l + H + C) A . (3.9) 

Hence, in (13. 5 p 

E\y{s + 9)- < [2^~'/5 + 3^-i2^-^iri(l + H{p)){l + H + C)A 

+ 2P-'E\yiv)-yik,A)\P. 



In what follows, we divide the following five cases to estimate the second term on the right- 
hand side of (13101) . 

Case 1: A;^ > and < - K/\ < A. By (EBD 

E\y{v)-y{K/\)Y 

= E\f{yk^A){v - KA) + g{yk^A){B{v) - B{KA)) + h{yk^A){N{v) - N{KA))\p 
< 3P-'E\fiy,^A)r{v - KAY + 3^-'E\giy,^A)\PE\Biv) - BiKA)]^ 
+ 3P-'E\h{yk^A)\PE\N{v) - N{KA)f. 

Then, in the same way as fl3.9p was done, we have for A G (0, 1) 

E\y{v) - yik,A)f < 3^^^2-2-^X2(1 + H{p)){l + H + C) A . 

Case 2: > and A < f — k^A < 2A. It then follows easily that 

E\y{v) - y{k,A)\P = E\y{v) - y{{kv + 1)A) + y{{kv + 1)A) - yik,A)\P 

< 2P-'E\yiv) - y{{kv + l)A)r + 2F-^E\y{{k, + 1)A) - y{k,A)\P. 

This, together with fl3.9p and Case 1, leads to 

E\y{v) - y{k^A)\P < 3^-^2^-^K2{l + H{p)){l + H + C) A . 

Case 3: ky = —1 and < f — kyA < A. In this case, —A < f < 0. We then have from 
(1231) that 

E\yiv)-y{kyA)\P<PA. 

Case 4- ky = —1 and A < v — k^A < 2A. In such case, < f < A. Case 1 and Case 2 can 
be used to estimate the term 

E\y{v) - y{KA)\^ < 2^-^E\y{v) - m\' + 2^-^E\m - y{{k,A)\^ 

< [2P-^f3 + 3P-^2'^-^K2{l + H{p)){l + H + C)]A. 

Case 5: ky < —2. In this case, t> < 0. So, by fl2.5p 

E\y{v)-y{k,A)f<2pA. 

Combining case 1 to case 5, we therefore complete the proof. □ 

The following Theorem will tell us the error of the pth moment between the true solution 
and numerical solution under global Lipschitz condition. 

Theorem 3.1. Under the conditions (12.31) and (12. 5p . forp > 2, 

E ( sup \x{t) - y{t)A < 6iL^e^^^^A, (3.11) 

\0<t<T J 

where 61,62 are constants which are independent of A. 
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Proof. It is easy to see from fl2.ip and fl2.8p that for any ti G [0, T] 



E ( sup \x{t) - y{t)\A < 3P-^E ( sup 

\0<t<ti / \0<t<ti 



f{xs-) - f{ys-)ds 



+ 3^-'e( snp [ g{xs-)-g{ys-)dB{s) \ 

\0<t<ti Jo J 

-^E ( sup / 

\0<t<ti Jo 



(3.12) 



+ 3^- 
:= h + l2 + h- 



h{x,-) - h{y,-)dN{s] 



In the sequel, we estimate these terms respectively. By the Holder inequality, (12.31) and 
Lemma [X! 



ti 



h<3^-'T^-' I E\f{x,)-fiy,)\^ds 
'o 



< &-'T^-' [ ' E\f{x,) - f{ys)\''ds + &-'T^-^ [ ' E\f{y,) - fiy^^ds 
Jo Jo 

p 

x{s + e) -y{s + 9)\'^dfx{9)] ' ds 



rti / pO 

Jo 



ti / i-O 





\y{s + e)-UO)rMO)] ds 



< QP-^TP-^ii e( sup \x{r) - y{r)A ds + Qp-^Tp-'l"^ [ ' [ E\y{s + 9)- ys{eWdfi{9)ds 

Jo Vo<r-<s / Jo J-T 

< &>-^TPL^-i A +QP-^Tp-^L'^ e( sup \x{r) - y{r)A ds. 

Jo \0<r<s J 



Now, the Burkholder-Davis-Gundy inequality |Fj Theorem 7.3, p40], (12. 3p and Lemma [3.21 
also give that, for some positive constant Cp, 



h<^'-'CpE \g{xs) - g{ys)\'ds 

< SP-'T'^Cp r E\g{xs) - g{ys)\'ds 
Jo 



< Qp-^T^a, 



T 2 



/ e( s\x 

Jo Vo<r 



sup \x{r) — y{r)\P ] ds + L 

0<r<s 

ti 



p-2 




E\y{s + 9) ~ yM^dmds 



J-T 



< QP-^-fT^CpL^ A +QP-^T'^CpL^ f ' E ( sup \x{r) - yir)A ds. 

Jo Vo<r<s / 



(3.13) 



In the same way as (I3.13P was done, together with the Burkholder-Davis-Gundy inequality 



pUl Theorem 48, pl93], we can deduce from fl2.3l) that, for some positive constant Cp, 



< QP-^E ( sup [ h{ 

\Q<t<ti Jo 



{x,-) - h{y,-)dN{s) 



+ XP sup 

0<t<ti 



h{xs) - h{ys)ds 



< QP-^CpT"^ X-^ + XPTP-^) I E\h{xs) - h{ys)\Pds 



< UP-^CpT'^X^ + XPTP~^)L2 e( sup \x{r) - y{r)A ds 

Jo Vo<r<s / 



+ 12P~\CpT'-^X^ + XpTP-^)L'^ J J E\y{s + 9) - y,{9)\Pdfi{e)ds 

<12P-^-f{CpT^X^ + XPTP)L^A+12P-\CpT^X^ + XPTP-^)L^ ['e[ sup \x{r) ~ y{r)\P ] ds. 

Jo Vo<r-<s 

Therefore 

E ( sup \x{t) - y{t)A < A +62L2 r e( sup \x{r) - y{r)\p\ ds, 

\0<t<ti J Jo \0<r<s J 

where Si = 6^-^7X5 (t^ +Cp + 2P-^CpX^ + X^T^^ and 62 = QT"^ (t§ +Cp + 2P-^CpX^ + 
A^Tty The desired assertion thus follows from the Gronwall inequality. 

us 



2-2 



ti rO 



□ 



Remark 3.1. The result of Theorem \3 . 1\ tells 



E ( sup \x{t)-y{t)A < SsLe^'^A, 

\0<t<T J 



(3.14) 



where 63, 64 are constants which are independent of A under the global Lipschitz condition 
(12.31) . This means that the order of the mean-square convergence is 1/2, while Eq. (13. lip 
tells us that the order of the pth moment convergence is \/p {p > 2). In other words, the 
lower moment has a better convergence rate for the SFDEs with jumps, whence it is best in 
practice to use the mean-square convergence. This is significantly different from the result 
on SFDEs without jumps. Letting h = in (12.11) . i.e. there is no jumps, we have already 
known that forp > 2 (see fT^ ) 



e( sup \xit) - yit)A <CiAP/\ 

\0<t<T J 



where Ci is a constant independent of A. This means that the order of the pth moment 
convergence is 1/2 for allp >2. Why is there a significant difference? Actually, it is due to 
the following fact: all moments of the Poisson increments N{{k + 1) A) — N{kA) have the 
same order of A (see {\JJ.8^ ) . while the moments of increments AB^ = B{{k -\-l)A) — B{kA) 



have different orders, namely E\ A Bf 



\2n 



0(A") and E\ A By 



\2n+l 



0. 
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4 Rate of convergence under local Lipschitz condition 



In this section, we shall discuss the rate of convergence of EM numerical solutions to 
(12.11) under the following local Lipschitz condition. 

(H3) (Local Lipschitz condition) For each integer j > 1, there is a left-continuous nonde- 
creasing function /ij : [— r, 0] i?+ such that 

- /wr V \9iv) - gm' V \hi^) - < f \m - m?dH{e), (4.1) 



for those V e D{[~t, 0]; i?") with \\^\\ V ||^|| < j. 

(H4) (Linear growth condition) Assume that there is a constant /i > such that, for ip G 
^([-r,0];i?"), 

\f{v)\'y\g{v)\'y\h{^)\'<h{i + \\^r). (4.2) 

Remark 4.1. Under the conditions (14. ip and (14.21) . for any initial data i e D^^^{ [-r, 0] ; i?") , 
(12.11) admits a unique solution x(t),t G [0,T] 6?/ using the standard truncation procedure 
(see /<§!, Theorem 3.4, p56] ). Moreover, (14. ip implies for those ip^ip E Z)([— r, 0]; i?") with 
IIV^II Vll^ll <3 

i/(^) - fm' V i^7(^) - gm' V iM^) - km' < h\w - n\ (4.3) 

where Lj = /ij(0) — fij{—r). 

Theorem 4.1. Let conditions (12.50 . (14. 1|) anc? (14.20 hold. If there exist positive constant a 
and 6 E (0, 1) such that the local Lipschitz constant obeys 

L]+'<a\ogj, (4.4) 

then 

E ( sup \x{t) - y{t)\A = O(A^), (4.5) 

\0<t<T J 

where e G (0, e) is an arbitrarily fixed small positive number. 

Proof. Let j > 1 be an integer, and let Sj = {x G : |x| < j}. Define the projection 
TTj : i?" ^ Sj by 



J A| 


x\ 




\x 





x, 



where we set 7rj(0) = as usual. It is easy to see that for all x,y E R 

kj(a:) - vrj(?/)| < \x - y\. 
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Define the operator tt^- : t){[-T, 0]; i?") D{[-t, 0]; i?") by 

7r,(y,) = {vr,(y,(^^)):-r<0<O}. 

Clearly, 

||7r,-(y.)|| <j, Vy.G^([-r,0];i?^). 

Define the truncation functions fj : i)([-r, 0]; ^ c/^ : Z)([-r, 0]; ^ i?'^^"" and 
/ij- : i^lhr, 0];i?") ^ i?" by 

= /(TTjM), ^i(<^) = ^(7ri(^)), hj{i^) = h{7ij{^)), (4.6) 

respectively. Then, by (14.11) . for any ip,ip E -D([— r, 0]; -R"), 

\fM - fM)? V \gM - g.m^ V M^) - /i,(^)r 

< |/(7r,(^)) - /(7r,(^))p V |(7(7r,(^)) - g{Tx,m? V |M7r,(^)) - W)P 

</° |v5(e)-^Wprf/i,W. 

That is, fj, Qj and hj satisfy the global Lipschitz condition. For t G [0,T], let x^{t) be the 
solution to the following SFDE with jumps 

dx^it) = fj{x{_)dt + gj{xi^)dB{t) + hj{x{_)dN{t) 

with the initial data Xg = ^ and i/-^(t) be the corresponding continuous-time EM solution 
with the stepsize A. By Theorem 13. II for any sufficiently smalle G (0,£) 

E ( sup \x^{t) - y'{t)\^^'] < 5iL;+^/'e^2^r'^' A . 

\0<t<T J 

Furthermore, by (14. 4p (here we assume Lj > 1 without any loss of generality), 

E ( sup \x^it) - y'{t)\^^'] < e^-'^+^^^^r^^'A < ^ _ (4^3) 

Set 

x(T) = sup \x(t)\ and y(T) = sup \y(t)\. 

0<t<T 0<t<T 

For any integer j > 1, define stopping time 

r, = T A inf{t G [0, T] : V \\yi\\ > j}. 
It is easy to see that || < j for any < s < Tj. Then, combining (14.61) gives that for any 

< S < Tj 
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Similarly, 

9j {x^s- ) = 9j+i {xl^ ) = 9{xl^ ) , hj {xl^ ) = hj+i {x^^ ) = h{xl. ) . 
While on < t < r, 

x'{t)=m+ f fj{xi-)ds+ f g,{xl_)dB{s) + f h,{xi_)dN{s) 
Jo Jo Jo 

= m+ [ fj+i{xi.)ds+ [ g,+i{xl.)dB{s) + [ h,+i{xl.)dN{s) 
Jo Jo Jo 

= m+ [ f{xi.)ds+ [ gixl.)dBis)+ [ hixl_)dNis). 
Jo Jo Jo 

Consequently, we must have that 

x{t) =x^{t) =x^+\t) 
on < t < Tj. Likewise, we can also derive that 



for < t < Tj. These imply that tj is non-decreasing and, by Lemma [3.11 lim Tj = T a.s. 
Let Tq = and compute, for t G [0, T], 

oo 
oo 



oo 



< \x^{t) -y\t)\'^I[j-l<£(T)Vy(T)]- 



Therefore, by the Holder inequality 



e( sup \x{t)-y{t)\^ 

\0<t<T 



(4.9) 



^ E (^sup^ \x'{t) - y'{t)\'+''^y^' [P{j - 1 < x{T) V y(T))]^. 

On the other hand, for any g > 2, we obtain from Lemma [3. II 

P(,-l<*(T)V«T))<^mM^,|M (,,0) 
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with j >2. Substituting (gS]) and fICTID into (HSj), one has 



(r,\ I <l<^ E ^ — ^ 2c<{(5i +<S2) — 9E \ 2 
sup \x{t) - y{t)\^ j < \^ + 2—^{2H{q))—^^j ^- j A^Ti . (4.II) 



For any fixed e > letting q be sufficiently large for 

+ ^2) + 2(2 + e) 



> 



we see that the right-hand side of fl4.11l) is convergent, whence the desired assertion fl4.5p 
follows. □ 

Remark 4.2. Under the local Lipschitz condition, Mao [9] showed the strong convergence of 
the numerical solutions to SFDEs without jumps, and the rate of convergence was revealed 
under the global Lipschitz condition. In the present paper, under the local Lipschitz condition, 
we reveal the rate of convergence for the numerical solutions to SFDEs with jumps. The rate 
of convergence for jump processes (12. ip we revealed here is 1/(2 + e) (closed to 1/2) under 
the logarithm growth condition (14. 4p . This is different from the rate of convergence for the 
diffusion processes (without jumps) which was studied in fW^ . where it was shown that the 
rate of convergence is still 1/2 under the logarithm growth condition. The reason for such a 
difference has already been pointed out in Remark 3.1. 



5 Appendix: an existence-and-uniqueness theorem 

To make our paper self-contained, in this section we shall discuss the existence and 
uniqueness of solutions to (12. ip under the assumption {HI). 

Theorem 5.1. Under the conditions (12. 2p . there exists a unique solution x(t),t G [0,T], to 
for any initial data ^ G D^^^{\-t, 0]; i?"). 

Proof. Since our proof is an application of the proof for the case without jumps in [8|, 
Theorem 2.2, pl50], we here give only a sketch for the proof of jump case. 

Uniqueness. Let x{t) and x{t) be two solutions to (12.11) on [0,T]. Noting from (12. ip that 

x{t)-x{t)= [ [f{xs-)-f{xs-)]ds+ [ [g{xs-)-g{xs-)]dB{s)+ [ [h{x,-) - h{x,-)]dN{s) 
Jo Jo Jo 

and N{t) = N{t) — At is a martingale measure for t e [0,T], along with (12. 3p we have 
E ( sup \x{t) - x{t)\A < 3L(T + 4 + 8A + 2X^T) [ E ( sup \x{r) - x{r)\A ds. 

\0<t<T J Jo \0<r<s J 

By the Gronwall inequality 

E I sup \x{t) -x{t)\A = 0, 
Vo<t<r / 
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which imphes that x{t) = x{t) for t G [0,T] almost surely. The uniqueness has been proved. 

Existence. Define Xq = ^ and x^{t) = ^{0) for < t < T. For each = 1, 2, ■ ■ ■ , set = ^ 
and define, by the Picard iterations, 

x-{t)=m+ ff{x:i')ds+ f g{x^z')dB{s) + f h{x^^-')dN{s) (5.1) 
Jo Jo Jo 

for t e [0,T]. It also follows from (15. ip that for any integer k>l 

E ( sup \x''{t)A <ci + C2 [ e( sup |x""i(r)p) ds, 

\0<t<T / Jo Vo<r<s / 

where Ci = 4[EUf+K{T^+AT) + {8X+2X'^T)T]+C2TEUf and Ca = 4K[T+4+8X + 2X^T]. 
This further implies that 

max E I sup |x"(t)|^ ) < ci + C2 / max E I sup \x"'~^{s)f ) dt. 

l<n<k \0<t<T J Jo \0<s<t J 

Observing 

max E I sup |a;"~"'^(s)p ) 

l<n<k \0<s<t J 

= max{E\^{0)\^,E ( sup \x\s)\A ,E ( sup \x''-\s)A\ 

I \0<s<t J \0<s<t J J 

< max 1^11^11', sup \x\s)A ,E ( sup \x''-\s)A , E ( sup \x\s)A\ 

I \0<s<i / \0<s<t J \0<s<t J ) 

< E\\if + max E ( sup |x"(s) 

l<n<k \n<s<t 



we hence deduce that 

max E ( sup \x''{t)A < Ci + C2TE\\^f + C2 I max E 

l<n<k \0<t<T J Jo 

The Gronwall inequality implies 

max E ( sup \x''{t)\A < {ci + C2TEUf)e'^^ . 

l<n<k \0<t<T J 

Since k is arbitrary, we must have for n > 1 

E ( sup \x^it)A < ic, + C2TEUf)e''^. 

\0<t<T ) 

Next, by dSlD 

^fsup |x^(t) -x°(t)|2^ < 3ir(T + 4 + 8A + 2A2r) /" {\^E\xlf)ds 

\0<t<T J Jo 

< 3KT{T + 4 + 8A + 2A2T)(1 + E^f) := C. 
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(5.2) 



We now claim that for n > 

sup - < , , < t < T, (5.3) 

o<s<t J n\ 

where M = 3K{T + 4 + 8A + 2A^T). We shall show this by induction. In view of fl^ we 
see that (15. 3p holds whenever n = 0. Now, assume that (15.31) holds for some n > 0. Then, 

E ( sup _ x''+\s)\A <M [ - x'^fds 

\0<s<t J Jq 



<m[ E f sup |x"+^(r) - x"(r)|2^ ds 

Jo \0<r<s J 

< M / : — ds 



n\ {n + l)\ 

Following the proof of [H Theorem 3.1, p55], we can show that for almost all a; G Q there 
exists a positive integer hq = uq^lj) such that 

sup \x"'~^'^{s) — a;"(s)| < — whenever n > uq^uj). (5.4) 

0<s<T 2" 

This implies {x"'{-)}n>i is a Cauchy sequence under sup | ■ |. However, since our space 
D{[0,T]; R"') is not a complete space under sup | ■ |, we do not know whether {x"'{-)}n>i 
has a limit in D{[0,T]; R^). In order for D([0, T]; i?") is complete, we need to define the fol- 
lowing metric (see [1, Chapter 3]). Let A denote the class of strictly increasing, continuous 
mapping of [0, T] onto itself and 



A* = <j A G A : sup 



< e 



define 



d{iX) = inf{e > : 3A G A* such that sup |^(t) - C(A(t))| < e}. 

te[o,T] 

d{-, ■) is called the Skorohod metric, by pj, Theorem 14.2, pll5] we know that -D([0, T]; i?") is 
complete in the metric d. Taking \{t) = t, we can see {x"'{-)}n>i is a Cauchy sequence under 
d. Therefore there exists unique xit),t G [0,T] G i5([0, T]; i?") such that d{x"'{-),x{-)) — > 
as n ^ oo. Taking the limit in (15.11) . we then can show that x{t) is the solution of (12.11) . □ 
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